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Abstract 

The  paper  will  provide  representative  simulations  of  particle  transport  around  a  vehicle  in  order  to  inves¬ 
tigate  some  of  the  issues  related  to  the  accurate  prediction  of  emission  and  transport  of  particles  induced 
by  a  moving  vehicle  with  a  transverse  blowing  wind.  Special  treatments  in  boundary  conditions  and  wall 
law  function  are  discussed  and  applied  to  maintain  the  shape  of  atmospheric  boundary  layer  wind  velocity 
profile.  For  the  vehicle,  we  adopt  the  geometry  of  a  Nissan  Pathfinder  SUV  to  study  the  effects  of  vehicle 
emission  and  transport  around  a  moving  vehicle.  We  perform  a  set  of  simulations  to  better  understand  the 
modeling  requirements  for  dust  emissions  including  a  sensitivity  study  to  determine  the  modeling  parame¬ 
ters  that  are  most  important  for  accurate  modeling  of  dust  generation  and  transport.  In  particular,  we  study 
the  effects  of  location,  size  distribution,  and  initial  velocity  distributions  of  the  modeled  dust  emissions  on 
predicted  downwind  atmospheric  dust  distributions. 

1  INTRODUCTION 

The  transport  of  dust  in  the  air  generated  by  vehicles  trav¬ 
eling  through  a  dusty  terrain  can  have  significant  impact 
on  many  aspects  of  vehicle  performance.  Operationally, 
dust  can  be  ingested  by  the  engine  and  power  train,  in¬ 
creasing  maintenance  expense  and  reducing  vehicle  life¬ 
time.  In  convoy  types  of  operations,  dust  emissions  can 
adversely  effect  visibility,  particularly  if  dust  accumulates 
on  vehicle  wind  shields  and  mirrors.  In  commercial  appli¬ 
cations,  the  transport  of  particulates  around  a  vehicle  is 
of  particular  interest  in  determining  vehicle  soiling. 

When  a  vehicle  travels  on  an  unpaved  road,  the  forces 
of  the  rolling  wheels  impacted  on  the  road  surface  crush 
large  particles  into  smaller  pieces  and  eject  them  into  the 
air.  The  vertical  pressure  of  the  wheel  on  the  ground 


causes  the  pulverization  of  surface  material.  Horizontal 
friction  which  sustains  the  vehicle  velocity  and  accelera¬ 
tion  of  the  vehicle  further  brings  the  particles  on  the  sur¬ 
face  of  the  tire.  The  slippage  between  the  tire  and  the 
road  lifts  particles  due  to  the  shear  force,  and  ejects  them 
into  the  air  due  to  centrifugal  force.  Saltation  is  another 
dust  generation  mechanism,  which  involves  the  bouncing 
of  bigger  particles  on  the  road,  and  thus  kicking  smaller 
particles  into  the  air.  Fine  particles  can  also  entrained 
in  the  turbulent  air  behind  the  vehicle.  Generally,  bigger 
particles  fall  back  to  the  ground  surface  more  rapidly  while 
smaller  particles  transfer  with  the  air  current.  Some  dis¬ 
cussions  on  dust  generation  mechanisms  can  be  found 
in  the  literature  [1,  2],  There  are  many  factors  affect  the 
generation  of  the  dust,  such  as  dust  particle  properties 
(size,  density,  stickiness,  etc.),  environment  conditions 
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(wind  speed, etc.),  soil  conditions  (silt  content,  moisture 
level,  etc)  and  vehicle  features  (vehicle  velocity,  weight, 
dimension,  number  of  the  wheels,  etc.).  There  have  been 
many  research  efforts  in  investigation  of  the  relationship 
between  the  magnitude  of  dust  emission  and  the  above 
factors  [3,  4,  5,  6],  However,  there  is  still  no  complete 
model  that  can  accurately  characterize  and  quantify  this 
relationship  due  to  the  complexity  and  variability  of  the 
problem. 

Most  experimental  work  on  dust  emission  induced  by 
vehicles  was  targeted  on  environmental  protection  (pol¬ 
lution  control)  applications  [6,  5,  7,  8,  9,  10],  thus  there 
is  a  lack  of  experimental  data  on  dust  concentration  in 
the  close  vicinity  of  vehicles.  Numerical  work  performed 
by  a  research  group  in  George  Mason  University  [2,  11] 
focused  on  a  real-time  simulation  of  dust  behavior  gen¬ 
erated  by  moving  vehicles.  Since  their  purpose  was  to 
simulate  the  dust  behavior  in  virtual  environment  for  in¬ 
teractive  graphic  applications,  and  the  computational  time 
had  to  be  fast  enough  to  realize  the  true  environment,  the 
numerical  models  were  over-simplified. 

In  this  study,  particle  transport  around  a  representa¬ 
tive  vehicle  (Nissan  Pathfound  SUV)  is  simulated  using 
Loci/CHEM  code  and  its  Lagrangian  particle  model.  A 
fully  developed  transverse  wind  profile  is  imposed.  In  the 
following  sections,  development  of  atmospheric  boundary 
layer  (ABL)  is  discussed  with  details  on  boundary  condi¬ 
tion  treatment  and  turbulence  model.  Then,  numerical 
simulations  on  particle  transport  in  the  air  induced  by  a 
moving  vehicle  will  be  presented,  with  the  discussion  of 
Lagrangian  particle  approach  and  preliminary  results  on 
sensitivity  of  the  dust  distribution. 


2  DEVELOPMENT  OF  ATMO¬ 
SPHERIC  BOUNDARY  LAYER 


function  which  is  based  on  boundary-layer  theory  for  the 
computational  cell  immediately  adjacent  to  the  wall.  We 
will  start  with  describing  the  wall  function  with  the  consid¬ 
eration  of  the  roughness  height.  The  turbulence  model 
adopted  in  this  study,  with  the  inclusion  of  wall  roughness 
treatment  will  be  presented.  To  evaluate  the  effective¬ 
ness  of  the  boundary  condition  treatment,  we  will  perform 
a  simulation  of  flow  over  a  flat  surface  with  one  kilometer 
in  length  to  determine  the  model’s  ability  to  preserve  the 
natural  ABL. 


2.1  Wall  Function  Boundary  Condition 

In  order  to  save  computational  expense,  a  wall  law  func¬ 
tion  is  often  prescribed  at  the  wall  boundary.  The  typi¬ 
cal  boundary  layer  velocity  profile  consists  of  a  thin  sub¬ 
layer  at  the  bottom  (y+  <  10),  a  logarithmic  profile  for  the 
30  <  y+  <  1000,  and  buffer  layer  in  between.  y+  is  the 
dimensionless  distance  to  the  wall  defined  as 


where  u*  and  v  are  friction  velocity  and  dynamic  fluid  vis¬ 
cosity  respectively.  There  are  several  empirical  relation¬ 
ships  available  for  the  velocity  in  a  boundary  layer  that 
can  be  used  to  develop  wall  function  boundary  conditions. 
The  formulation  of  Spalding  [13]  is  adopted  in  this  study, 
and  it  is  a  unified  form  valid  for  the  log  layer  and  sublayer 
as  well  as  the  transition  region.  This  form  is  given  by 


y+  =u+  +  e-KB[eKU 


—  1  —  Kir 


(kw+)2 

2 


(kw+)3 
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(2) 


in  which  k  and  B  are  dimensionless  constants  in  the  log 
layer  profile  (k  =  0.41),  u+  is  defined  as 


Since  dust  particles  are  transported  and  diffused  in  the 
air,  it  is  necessary  to  accurately  describe  the  atmospheric 
boundary  layer.  Because  the  simulation  domain  is  nor¬ 
mally  much  smaller  compared  to  the  scale  of  the  terrain, 
it  is  not  possible  to  obtain  a  fully  developed  ABL  by  ap¬ 
plying  a  uniform  flow  at  the  domain  boundaries.  Thus, 
a  prescribed  ABL  velocity  profile  is  normally  set  at  the 
inflow,  and  it  is  important  that  this  profile  maintains  its 
shape  along  the  fetch.  Richard  and  Hoxey  [12]  have  pro¬ 
posed  a  set  of  boundary  conditions  to  ensure  a  homoge¬ 
neous  boundary  layer.  The  standard  logarithmic  velocity 
profile  with  the  consideration  of  surface  roughness  length 
is  imposed  on  the  inflow  boundary.  At  the  top  boundary, 
a  constant  shear  stress  is  applied  to  compensate  the  en¬ 
ergy  lost  due  to  the  shearing  on  the  ground.  In  addition, 
the  surface  roughness  has  to  be  incorporated  either  in 
the  boundary  condition  of  viscous  wall,  or  through  a  wall 


u 


+ 


u 

U * 


(3) 


For  smooth  wall,  the  constant  B  is  approximately  set  to 
5.2.  For  surfaces  with  average  roughness  height  of  ks, 
the  above  wall  law  still  holds.  However,  the  additive  con¬ 
stant,  B  is  a  function  of  ks  in  the  rough  wall  cases.  Cebeci 
and  Bradshaw  [14]  suggested  the  following  correlation  of 
Nikuradse’s  data  [15]  to  evaluate  the  B: 


(—3.3  +  ^fnK+ )  sin[0.4258(fnK+— 
0.8111)] +  5.2, 

(-3.3+  \lnn+)  +5.2,  k+>90 


2.25  <  k+  >  90  (4) 


This  correlation  is  implemented  to  compute  the  constant 
B  if  wall  roughness  is  considered. 
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2.2  Baseline  Turbulence  Model  (BSL) 

There  are  several  turbulence  models  implemented  in  the 
CHEM  solver.  Since  the  wall  roughness  treatment  is  only 
implemented  in  the  Menter’s  Baseline  model  (BSL)  and 
Wilcox  (1998)  [16]  models,  and  BSL  model  is  an  modifi¬ 
cation  of  Wilcox  (1998)  model,  we  will  only  describe  the 
BSL  model  in  the  following. 

It  is  well  known  that  two-equation  eddy-viscosity  low- 
Reynolds-number  turbulence  models  are  among  the 
most  widely  used  models  for  engineering  applications  to¬ 
day,  and  the  k  -  e  model  with  damping  functions  near 
the  wall  is  the  most  popular.  However,  the  k  -  e  model 
often  suffers  from  numerical  stability  problems  due  to  dis¬ 
parate  turbulent  time  scales.  Another  well-known  two- 
equation  turbulence  model  is  the  k  -  ui  model,  developed 
by  Wilcox[17],  It  has  the  advantage  that  it  does  not  re¬ 
quire  damping  functions  in  the  viscous  sublayer  and  that 
the  equations  are  less  stiff  near  the  wall,  therefore  it  is 
superior  to  k  -  e  model  with  regard  to  numerical  stabil¬ 
ity.  However,  when  applied  to  the  free  shear  layers,  it  is 
found  that  there  is  a  strong  dependency  of  the  results  on 
the  free  stream  value  of  w[18,  19].  Menter  created  a  new 
model,  called  baseline  (BSL)  model,  by  blending  the  k  -  e 
and  k  -  oj  models[20].  It  utilizes  the  k  -  ui  model  in  the 
wall  region  and  gradually  switches  to  the  k-e  model  away 
from  the  wall.  To  achieve  this,  the  k-e  model  is  first  trans¬ 
formed  into  a  k  -  oj  formulation,  and  an  additional  cross 
diffusion  term  is  added  (another  diffusion  term  associated 
with  turbulent  kinetic  energy  is  neglected  in  the  formula¬ 
tion  under  certain  assumptions[16]).  The  original  k  -  ui 
equations  are  then  multiplied  by  a  blending  function  FbsU 
the  transformed  k-e  equations  are  multiplied  by  (1  -Fbsi), 
and  then  both  are  added  together.  The  blending  function 
Fbsi  is  designed  so  that  it  is  unity  at  the  wall,  and  grad¬ 
ually  approaches  zero  away  from  the  wall.  Note  that  the 
k  -  lo  model  can  be  easily  obtained  by  setting  Fbsi  =  1 
identically. 

The  defining  equations  for  the  BSL  model  are  written 
as: 

Kinematic  Eddy  Viscosity: 


vt  =  k/ui ,  (5) 

where  vt  is  the  turbulent  viscosity,  k  is  the  turbulence  ki¬ 
netic  energy,  and  ui  represents  the  dissipation  per  unit 
turbulence  kinetic  energy. 

Turbulent  Stress  Tensor: 

/  ( dui  duj  \  2  _ 

i  =  1, 2, 3,  j  =  1,2,3, 


Turbulent  Kinetic  Energy  Equation: 


Dk  •  dui 

P^t=Tijd^ 


(3*  puk  +  -g—  [fa  +  ptCTk) 


dk  - 
dxjF 
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Turbulent  Dissipation  Equation: 

Dui 


Dt 


7  '  dui  2 

-T«aiT's'" 


d  - 

dxj  . 


dui 


,  .  1  dk  dui 

-(1  -  Fbsl)p<Ju: 2-faTfaE- 

UJ  1/tLj  L/J/j 


Closure  Coefficients: 

All  the  constants  <j>  of  the  model  are  computed  by  blend¬ 
ing  the  appropriate  k-ui  and  k-e  constants,  as  follows: 


<t>  =  Fbsifa  +  (1  —  Fbsi)<j>2,  (7) 

where  the  constants  fa  (k  -  ui)  are  used  are  provided  by 
Menter[20] 

The  blending  function  Fbsi  is  defined  as  follows: 


Fbsi  =  tanh(argbsl ), 


(8) 


where 


argbsi  =  mm 


max 


(: 


Vk  500r/\  Apau2k  1 

f^)'CDkuyfa'  [  ’ 


V0.09  uoy  y 2 

and  y  is  the  distance  to  the  closest  point  away  from  the 
wall  surface.  In  the  above,  CDkui  is  defined  as: 

1  dk  dui 


/  1  OK  OUI  _ on 

CDku;  =  max(2pcr^2  —  ^ — tt — ,10 
V  ui  dxj  dxj 


(10) 


The  boundary  conditions  for  k  and  w  at  a  solid  wall  are: 


k  =  0, 


lo  =  10 


ft  (A yi)2  ’ 


(11) 


where  Ayi  is  the  distance  from  the  first  cell  center  to  the 
solid  wall.  For  rough  walls,  the  user-specified  equivalent 
roughness  height  ( ks )  is  used  to  compute  the  solid  wall 
value  of  to  from: 

u2 

oi  =  ^Sr  (12) 

v 

where 


SR  = 


k+  <  5 


100 

kf 


K>  s 


(13) 


A  surface  is  considered  to  be  hydraulically  smooth1  when 
k+  =  uTks/v  <  5.  A  slightly  rough  wall  boundary  condi¬ 
tion  for  ui  can  be  derived  by  combining  Eqs.  (1 2)  and  (1 3) 
and  to  obtain[21  ] : 


40000^ 


(14) 


in  which  ks  must  be  chosen  to  ensure  that  <  5. 

1 A  hydraulically  smooth  surface  is  one  in  which  the  roughness  height 
is  smaller  than  the  laminar  sublayer. 
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2.3  Inflow  and  Top  Boundary  Conditions 

To  test  our  ability  to  simulate  the  ABL  we  adopted 
Richards  and  Hoxey’s  [12]  inflow  and  top  boundary  con¬ 
ditions  in  our  study  since  some  aspects  of  their  work  have 
become  a  standard  in  computational  wind  engineering  in 
order  to  sustain  the  ABL  velocity  profile.  Based  on  the 
k  —  e  turbulence  model,  Richards  and  Hoxey  suggested 
the  profiles  for  velocity,  turbulent  kinetic  energy  and  tur¬ 
bulent  dissipation  which  satisfy  the  the  k  -  e  model.  The 
profile  functions  are  written  as 


(15) 

«  V  yo  J 

9. 

II 

-Sc: 

(16) 

< 

K(y  +  yo) 

(17) 

where  C M  is  a  constant  in  k  -  e  model,  equal  to  0.09. 
Note  that  y0  is  the  surface  roughness  length  which  is 
different  from  surface  roughness  height,  ks.  There  are 
numerous  work  on  the  relationship  between  ks  and  y0 
[22](Hargreaves).  If  the  eq.(15)  is  rewritten  in  the  form 
of 

u+  =  —  ln(y+)  +  B  (18) 

K 

and  is  related  to  the  formulation  of  the  constant  consid¬ 
ering  the  roughness,  eq.(4),  we  have  ks  «  30 y0  and 
B  =  -14.9  in  the  current  study. 

The  profiles  of  velocity,  turbulent  kinetic  energy  and  tur¬ 
bulent  dissipation  represented  through  eq.(1 5)  to  (1 7)  are 
imposed  in  the  inflow.  At  the  top  of  computational  domain, 
a  constant  shear  stress  of  pul  is  applied,  as  proposed  by 
Richard  and  Hoxey,  to  maintain  the  energy  balance. 

2.4  Results  on  Boundary  Layer  on  a  Flat 
Plate 

Air  with  far  field  temperature  of  293AT  and  pressure  of 
1  atm  flowing  over  a  flat  plate  is  carried  out  to  test  the 
homogeneity  of  velocity  field  in  the  flow  direction.  In 
our  test  case,  we  use  the  same  prescribed  inflow  veloc¬ 
ity  profile  as  in  Richards  and  Hoxey’s  [12]  work.  Thus, 
the  friction  velocity  (m*)  and  the  roughness  length  (y0)  in 
the  eq.(15)  are  set  to  be  0.6 m/s  and  0.01m  respectively, 
which  were  obtained  by  fitting  the  measured  velocity  data. 
The  ground  roughness  height,  ks  is  set  to  0.3?n  according 
to  the  relationship  ks  «  30 y0,  which  leads  to  =  12, 000. 
Roughness  length  (y0)  of  0.01m  corresponds  to  the  ter¬ 
rain  condition  of  a  lawn.  For  a  level  desert,  y0  is  on  the 
order  of  0.001m.  Wilcox  (1998)  turbulence  model  [16]  is 
used  in  this  simulation  since  we  found  that  Wilcox  (1998) 


model  converges  better  that  BSL  model  for  roughness 
height  larger  than  0.1m.  The  wall  law  boundary  condition 
is  applied  on  the  ground  with  y+  equal  to  150  at  the  first 
grid  point  away  from  the  wall.  The  computational  domain 
is  1  km  x  0.5 km  in  flow  and  vertical  direction  respectively. 

The  development  of  boundary  layer  along  the  plate  is 
shown  in  figure  1 .  At  first  sight,  it  seems  that  the  velocity 
profile  maintains  itself  well  along  the  domain.  However, 
by  zooming  into  the  bottom  of  10m  of  boundary  layer,  the 
discrepancy  between  the  original  inflow  profile  and  down¬ 
stream  profile  becomes  apparent,  as  shown  in  figure  2. 
Figure  2  also  includes  the  results  without  the  roughness 
treatment  in  the  wall  function  and  shear  stress  applica¬ 
tion  at  the  top.  It  is  clear  that  the  velocity  profile  main¬ 
tains  a  much  closer  shape  to  the  prescribed  inflow  profile 
if  the  roughness  height  is  considered  in  the  wall  function 
and  the  shear  stress  is  imposed  at  the  top.  There  are  a 
few  possible  reasons  that  can  be  contributed  to  the  non¬ 
sustainability  of  the  ABL  even  with  Richards  and  Hoxey 
treatments.  First,  Richards  and  Hoxey  model  was  based 
on  k  —  e  model,  while  k  —  u  model  is  adopted  to  compute 
the  turbulence  in  this  study,  which  could  cause  some  non¬ 
consistency  regarding  the  model  parameters.  Second, 
there  is  difference  between  the  prescribed  friction  veloc¬ 
ity  and  actual  friction  velocity  calculated  from  the  code,  as 
depicted  in  figure  3.  Lower  friction  velocity  results  in  the 
larger  flow  speed  at  the  bottom  of  boundary  layer,  which 
is  seen  in  the  figure  2.  The  parameter,  B,  in  the  wall  func¬ 
tion  is  based  on  the  prescribed  friction  velocity,  and  re¬ 
mained  constant  during  the  development  of  the  boundary 
layer.  The  iterative  calculation  between  computed  friction 
velocity  and  the  corresponding  parameter,  B  in  the  wall 
function  is  proposed  to  eliminate  the  difference  between 
the  targeted  and  real  friction  velocity.  Maintaining  ABL 
profile  is  not  a  trivial  problem,  and  more  effort  will  be  car¬ 
ried  on  in  the  future  research.  In  the  current  application 
where  computational  domain  is  much  smaller  than  1  km, 
the  sustainability  of  the  ABL  is  satisfactory. 

3  PARTICLE  TRANSPORT  IN  THE 
AIR  INDUCED  BY  A  MOVING  SUV 

In  this  study,  we  have  performed  a  number  of  represen¬ 
tative  simulations  of  particle  transport  around  a  repre¬ 
sentative  vehicle  in  order  to  investigate  some  of  the  is¬ 
sues  related  to  predicting  emissions  and  transport  of  par¬ 
ticles  from  a  vehicle  traveling  with  a  transverse  wind.  For 
these  simulations  we  utilize  the  Loci/CHEM  code  and 
its  Lagrangian  particle  transport  model.  For  the  vehicle 
we  adapted  the  geometry  of  a  Nissan  Pathfinder  SUV 
that  was  previously  used  in  blast  studies  in  the  SIMBRS 
WD19  task.  A  viscous  mesh  was  generated  on  this  SUV 
using  the  SolidMesh/AFLR3  software  developed  at  Mis- 
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U  velocity  (m/s) 

Figure  1 :  X-component  velocity  distribution  in  the  vertical 
direction  along  the  flat  plate 


sissippi  State.  The  final  mesh  was  about  15  million  cells. 
Cases  were  run  on  48  processors  of  our  cluster  with  turn 
around  times  on  the  order  of  1-2  hours.  The  Lagrangian 
particle  model,  provided  by  a  scalable  parallel  implemen¬ 
tation  using  the  Loci  framework,  is  described  in  the  fol¬ 
lowing  section. 

3.1  Lagrangian  Particle  Approach 

The  governing  equations  Eqns.  (19-20)  for  the  particle 
movements  were  derived  using  the  Basset-Boussinesq- 
Oseen  (BBO)  assumption  that  the  density  of  the  particle 
is  much  larger  than  that  of  the  fluid  and  particle  size  is 
small  compared  to  turbulence  integral  length  scale,  and 
that  the  effect  of  shear  on  particle  motion  is  negligible: 


dxp 

dt 


(19) 


=  Dp(u  -  Up)  +  Fv  (20) 

where  up  is  the  particle  velocity,  and  u  is  the  gas-phase 
velocity  interpolated  to  the  particle  location  xp.  Fv  is 
the  body  force  which  includes  the  gravity.  Only  momen¬ 
tum  coupling  is  considered  in  the  current  study.  The 
drag  force  on  a  solid  particle  is  modeled  using  a  drag- 
coefficient, 


3  p  |  u  up  | 

=  - 3 - 

4  pp  d 


(21) 


Figure  2:  X-component  velocity  distribution  in  the  vertical 
direction  along  the  flat  plate  in  the  bottom  10m  of  bound¬ 
ary  layer 


where  p  and  pv  are  the  density  of  fluid  and  the  particle 
respectively,  dp  is  the  particle  diameter  and  Cd  is  obtained 
from  the  nonlinear  correlation  [23]: 


Figure  3:  Comparison  of  prescribed  friction  velocity  and 
the  computed  one 


Cd  =  j^-(l  +  aRebp)  (22) 

where  Rep  =  dp\u  -  up\/ p  is  the  particle  Reynolds  num¬ 
ber.  The  constants  a  =  0.15,  b  =  0.687  are  used  to  get 
the  drag-coefficient. 

The  Lagrangian  particle  implementation  uses  an  im¬ 
plicit  backward  differentiation  scheme  to  support  stiff 
source  terms.  The  scheme  is  expressed  in  a  generic  form 
as: 

s 

XP+1  +  E  =  A//3<+1  (23) 

i= 1 


<+1  +  E^uP+1_i  =  A t/3Dp(un+1  -  u%+1) 

i=  1 

+  AtpF™+1 
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where  s  is  the  total  number  of  the  multi-step  scheme, 
At  is  the  time  step,  and  coefficients  for  the  second  or¬ 
der  scheme  used  in  current  implementations  are:  a\  = 

-|,a2  =  and  (3  =  -§. 

The  injection  of  Lagrangian  particles  from  the  injection 
surface  is  needed  to  be  handled  in  a  more  realistic  way  in 
the  Lagrangian  approach.  The  size  and  location  of  each 
injected  particle  is  random.  Some  possible  particle  size 
distribution  functions  are  the  logarithmic  normal  distribu¬ 
tion  formulated  as 

Pr(d)  =  [1/ a  d]exp{— (\  /  2<r2)[ln{d  /  dmed)}2} 

given  the  median  particle  diameter  dmed  and  a  standard 
deviation  a,  and  the  skewed  logarithmic  distribution: 

Pr{d )  =  [1  -  (d/dmax)2][l  -  {dmin/d)2} 
x  exp{— (l/2a2)[ln(d)  —  /3]2} 

where  j3  is  computed  as 

2  a2 

/3  =  ln{dpeak )  +  — - — - -T— 

y^max  /  dpeak  ) 

given  peak, minimum,  and  maximum  particle  diameter 
dpeak,  dmin,  dmax ,  respectively,  and  a  standard  deviation 

a. 

Since  the  flow  field  information  is  needed  at  the  particle 
location  during  the  Lagrangian  particle  equation  integra¬ 
tion  process,  it  is  necessary  to  find  in  which  control  vol¬ 
ume  the  particle  resides.  Thus,  the  particle  search  and 
location  schemes  and  interpolation  schemes  are  imple¬ 
mented. 

The  purpose  of  the  algorithm  is  to  determine  if  a  par¬ 
ticle  exists  in  a  given  control  volume  with  arbitrary  shape 
[24,  25,  26,  27,  28].  For  the  CHEM  particle  simulation  we 
use  the  projection  approach  as  utilized  by  Apte  et  a/[28j. 
To  enhance  robustness,  we  mark  particles  that  have  an 
extended  path  length  and  track  the  cells  through  which 
they  pass.  Within  any  single  time  step,  we  do  not  allow  a 
particle  to  pass  through  the  same  cell  twice  while  search¬ 
ing  for  its  containing  cell.  This  produces  a  robust  geomet¬ 
ric  searching  algorithm  for  locating  a  particle  within  the 
cell  of  arbitrary  polyhedral  unstructured  meshes. 

The  flow  properties  are  interpolated  using  the  limited 
cell  reconstructions  from  the  fluid  phase  simulations.  This 
provides  a  spatially  second  order  interpolation  method. 

3.2  Results  on  Particle  Transport 

For  the  simulation  of  particle  generation  and  transport 
we  consider  an  SUV  type  vehicle  traveling  with  13.5m/s 
(about  30mph)  uniform  forward  velocity  with  a  cross  wind 
of  6.75  m/s  (about  15mph).  The  wind  is  simulated  either 


as  a  uniform  flow  or  an  ABL  with  a  velocity  of  6.75m/s 
at  one  meter  height  using  a  surface  friction  velocity  of 
w*  =  0.6m/ s  and  roughness  length  of  y0  =  0.01m.  The 
rubber  part  of  the  tire  is  prescribed  a  rotational  velocity 
348.42  rpm  which  is  consistent  with  a  forward  velocity  of 
1 3.5m/s.  The  rotation  of  the  spokes  of  the  hub  is  not  con¬ 
sidered  in  this  model. 


Particles  are  emitted  using  a  tri-mode  particle  distribu¬ 
tion  with  particle  mean  diameters  of  2.5 pm,  10 pm  and 
100/i?7t  with  a  mass  distribution  of  1 0%,  80%,  and  1 0% 
respectively.  Each  size  group  is  emitted  with  a  log-normal 
distribution  of  its  mean  particle  diameter  with  a  variance 
of  a  =  0.25.  Particles  have  a  material  density  of  quartz 
which  is  2648  kg/m3  and  are  injected  from  the  tires  at 
a  mass  flux  of  0.1  kg/m2  with  a  speed  of  95%  of  the 
wheel  velocity  for  the  cases  where  tires  were  treated  as 
the  source  and  a  similar  mass  flux  from  a  small  patch 
under  the  tires  at  a  fixed  upward  velocity  of  4m/s  in  one 
case.  The  BSL  turbulence  model  was  employed  in  all  sim¬ 
ulations.  The  effects  of  turbulence  diffusion  on  the  parti¬ 
cles  is  not  considered  in  this  simulation  as  we  will  need 
to  implement  a  turbulence-particle  coupling  term.  We  will 
likely  use  a  first  order  continuous  Markov  model  to  model 
turbulent  velocity  fluctuations  in  future  simulations. 


Our  first  simulation  considers  a  prescribed  uniform  flow 
for  the  wind  with  particles  emitted  from  the  tire  surfaces. 
A  side,  back,  and  top  view  of  this  simulation  is  shown 
in  figures  4,  5,  and  6  with  the  vehicle  surface  colored 
by  computed  pressures.  Particles  are  colored  by  a  gray 
scale  indicating  particle  diameter  with  larger  particles  be¬ 
ing  darker  in  color.  These  view  show  that  the  particle 
model  qualitatively  behaves  as  expected,  with  particles 
being  entrained  in  the  wake  flow  of  the  vehicle  much  as 
one  would  expect.  We  then  compare  changing  the  loca¬ 
tion  of  particle  emissions  from  the  tires  to  the  ground  in  a 
small  region  around  the  tire.  A  front  view  of  the  particle 
distributions  for  the  uniform  flow  showing  the  tire  emission 
is  illustrated  in  figure  7  while  the  results  from  the  ground 
emission  is  shown  in  figure  8.  Although  the  ground  emis¬ 
sions  result  in  a  particulate  cloud  that  is  hangs  closer  to 
the  ground,  the  two  cases  are  surprisingly  similar  sug¬ 
gesting  that  there  is  only  modest  sensitivity  to  the  means 
in  which  particles  are  injected  into  the  simulation.  Finally, 
we  consider  the  effects  of  the  atmospheric  boundary  layer 
on  the  simulation  (this  time  with  emissions  from  the  tires) 
and  find  that  there  are  significant  differences  in  the  plume 
shape.  We  expect  that  these  differences  will  become 
more  significant  as  we  add  turbulent  effects  to  the  par¬ 
ticle  coupling  as  the  atmospheric  boundary  layer  will  sig¬ 
nificantly  increase  particle  scattering  behind  the  vehicle. 
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Figure  4:  Side  view  of  particle  distribution  for  tire  emission  using  uniform  wind  model.  Vehicle  colored  by  pressure. 


Figure  5:  Back  view  of  particle  distribution  for  tire  emission  using  uniform  wind  model.  Vehicle  colored  by  pressure. 
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Figure  6:  Top  view  of  particle  distribution  for  tire  emissions  using  uniform  wind  model.  Vehicle  colored  by  pressure. 


Figure  7:  Front  view  of  particle  distribution  for  tire  emissions  using  uniform  wind  model.  Vehicle  colored  by  pressure. 
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Figure  8:  Front  view  of  particle  distribution  for  ground  emissions  using  uniform  wind  model.  Vehicle  colored  by 
pressure. 


Figure  9:  Front  view  of  particle  distribution  for  tire  emissions  using  atmospheric  boundary  layer  wind  model.  Vehicle 
colored  by  pressure. 


Sensitivity  Studies  on  Particle  Emissions  and  Transport  Around  a  Moving  Vehicle 


UNCLASSIFIED 

Proceedings  of  the  2011  Ground  Vehicle  Systems  Engineering  and  Technology  Symposium(GVSETS) 


Disclaimer 

Reference  herein  to  any  specific  commercial  company, 
product,  process,  or  service  by  trade  name,  trademark, 
manufacturer,  or  otherwise,  does  not  necessarily  consti¬ 
tute  or  imply  its  endorsement,  recommendation,  or  favor¬ 
ing  by  the  United  States  Government  or  the  Department 
of  the  Army  (DoA).  The  opinions  of  the  authors  expressed 
herein  do  not  necessarily  state  or  reflect  those  of  the 
United  States  Government  or  the  DoA,  and  shall  not  be 
used  for  advertising  or  product  endorsement  purposes. 


References 

[1]  M.  Succarieh.  Control  of  dust  emissions  from  un¬ 
paved  roads.  Technical  report,  University  of  Alaska, 
May  1992.  Prepared  for  Alaska  Cooperative  Trans¬ 
portation  and  Public  Facilities  Research  Program. 

[2]  J.  Chen,  X.  Fu,  and  E.  Wegman.  Real-time  simu¬ 
lation  of  dust  behavior  generated  by  a  fast  traveling 
vehicle.  ACM  Transactions  on  Modeling  and  Com¬ 
puter  Simulation ,  9(2)  :81  -1 04,  1 999. 

[3]  C.  Cowherd.  Sandblaster  2  support  of  see-through 
technologies  for  particulate  brownout.  Technical  re¬ 
port,  2007.  Task  5  Final  Technique  Report,  MRI 
Project  No.  110565. 

[4]  US  Environmental  Protection  Agency,  Research  Tri¬ 
angle  Park,  NC.  AP-42,  Miscellaneous  Source 
Chapter  13,  fifth  ed.  edition,  2003. 

[5]  J.  Gillies,  V.  Etyemezian,  H.  Kuhns,  H.  Mmsmuller, 
J.  Engelbrecht,  J.  King,  S.  Uppapalli,  G.  Nikolic, 
J.  McAlpine,  D.  Zhu,  M.  Skiba,  and  D.  Gillette. 
Particulate  matter  emissions  factors  for  dust  from 
unique  military  activities.  Technical  report,  DRI, 
2010.  SERDP  Project  SI-1399. 

[6]  J.  Gillies,  V.  Etyemezian,  H.  Kuhns,  D.  Nikolic,  and 
D.  Gillette.  Effect  of  vehicle  characteristic  on  un¬ 
paved  road  dust  emissions.  Atmospheric  Environ¬ 
ment,  39:2341-2347,  2005. 

[7]  J.  Veranth,  E.  Pardyjak,  and  G.  Seshadri.  Vehicle¬ 
generated  fugitive  dust  transport:  analytic  models 
and  field  study.  Atmospheric  Environment,  2003. 

[8]  H.  Kuhns,  V.  Etyemezian,  and  P.  Shinbein.  Relat¬ 
ing  road  dust  emissions  surrogates  to  average  daily 
traffic  and  vehicle  speed  in  las  vegas,  nevada.  Tech¬ 
nical  report,  DRI  and  Clark  County  Comprehensive 
Planing,  not  know. 


[9]  V.  Etyemezian,  H.  Kuhns,  J.  Gilles,  M.  Green, 

M.  Pitchford,  and  J.  Watson.  Vehicle-based  road 
dust  emission  measurement:  l-methods  and  cal¬ 
ibration.  Atmospheric  Environment,  pages  4559- 
4571, 2003. 

[10]  V.  Etyemezian,  H.  Kuhns,  J.  Gilles,  J.  Chow,  K.  Hen¬ 
drickson,  M.  McGown,  and  M.  Pitchford.  Vehicle- 
based  road  dust  emission  measurement:  lii-effect  of 
speed,  traffic  volume,  location,  and  season  on  pmw 
road  dust  emissions  in  the  treasure  valley,  id.  Atm- 
sopheric  Environment,  pages  4583-4593,  2003. 

[11]  J.  Chen,  E.  Wegman,  and  J.  Wang.  Animation  of 
dust  behaviors  in  a  networked  virtual  environment. 
This  article  is  found  using  google.  No  information  on 
the  year  and  where  published. 

[12]  P.  Richards  and  R.  Hoxey.  Appropriate  boundary 
conditions  for  computational  wind  engineering  mod¬ 
els  using  the  k-e  turbulence  model.  Journal  of  Wind 
Engineering  and  Industrial  Aerodynamics,  1 993. 

[13]  D.  B.  Spalding.  A  single  formular  for  the  law  of  the 
wall.  Transactions  of  the  ASME,  Series  E:  Journal  of 
Applied  Mechanics,  5(1 7),  1 961 . 

[14]  T.  Cebeci  and  P.  Bradshaw.  Momentum  Transfer  in 
Boundary  Layers.  Hemisphere  Publishing,  1977. 

[15]  J.  Nikuradse.  Laws  of  flow  in  rough  pipes.  Technical 
Report  Tech.  Rep.  NACA  TM  1292,  NACA,  1950. 

[16]  D.  C.  Wilcox.  Turbulence  Modeling  for  CFD.  DCW 
Industries,  1998. 

[17]  D.  C.  Wilcox.  Reassessment  of  the  scale¬ 
determining  equation  for  advanced  turbulence  mod¬ 
els.  Al AA  Journal,  26(1 1):1 299-1 310,  1988. 

[1 8]  D.  C.  Wilcox.  A  half  century  review  of  the  k—u  model. 
Technical  report,  AIAA,  1991.  AIAA91-1784. 

[19]  F.  R.  Menter.  Influence  of  freestream  values  on 
k  -  uj  turbulence  model  predictions.  AIAA  Journal, 
30(6):1 657-1 659,  June  1992. 

[20]  F.  R.  Menter.  Two-equation  eddy-viscosity  turbu¬ 
lence  models  for  engineering  applications.  AIAA 
Journal,  32(8):1 598—1 605,  August  1994. 

[21]  D.  C.  Wilcox.  Formulation  of  the  k  -  oj  turbulence 
model  revisited.  AIAA  J.,  46(1 1):2823-2838,  2008. 

[22]  J.  Franke,  C.  Hirsch,  A.  G.  Jensen,  H.  W.  Schatz- 
mann,  P.  S.  Westbury,  S.  D.  Miles,  J.  A.  Wisse,  and 

N.  G.  Wright.  Recommendations  on  the  use  of  cfd 
in  wind  engineering.  In  J.P.A.J  van  Beeck,  editor, 
Proceedings  of  the  International  Conference  Urban 


Sensitivity  Studies  on  Particle  Emissions  and  Transport  Around  a  Moving  Vehicle 


UNCLASSIFIED 

Proceedings  of  the  2011  Ground  Vehicle  Systems  Engineering  and  Technology  Symposium(GVSETS) 


Wind  Engineering  and  Building  Aerodynamics,  von 
Karman  Institute,  2004. 

[23]  M.  Sommerfeld  C.  Crowe  and  Y.  Tsuji.  Multiphase 
flows  with  droplets  and  particles.  CRC  Press,  1998. 

[24]  R.  Lohner  and  J.  Ambrosiano.  A  vectorized  particle 
tracer  for  unstructured  grids.  Journal  of  Computa¬ 
tional  Physics,  91 :22-31 ,  1 990. 

[25]  T.  Westermann.  Localization  schemes  in  2d 
boundary-fitted  grids.  Journal  of  Computational 
Physics,  101:307-313,  1992. 

[26]  R.  Lohner.  Robust,  vectorized  search  algorithms  for 
interpolation  on  unstructured  grids.  Journal  of  Com¬ 
putational  Physics,  118:380-387,  1995. 

[27]  A.  Allievi  and  R.  Bermejo.  A  generalized  particle 
search-locate  algorithm  for  arbitrary  grids.  Journal 
of  Computational  Physics,  1 32:1 57-1 66,  1 997. 

[28]  P.  Moin  S.  V.  Apte,  K.  Mahesh  and  J.  C.  Oe- 
felein.  Large-eddy  simulation  of  swirling  particle¬ 
laden  flows  in  a  coaxial-jet  combustor,  international 
Journal  of  Multiphase  Flow,  29:1311-1331,  2003. 


Sensitivity  Studies  on  Particle  Emissions  and  Transport  Around  a  Moving  Vehicle 


